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Statement  of  problem  studied 

The  principal  goal  of  this  research  was  to  develop  a  system  to  experimentally  measure  the  three-dimensional  transient 
deformation  processes  that  occur  in  multi-layer,  functionally  graded  ceramic  (CMC’s)  armor  systems  undergoing  a  range  of  dynamic 
loading  conditions.  The  development  of  the  metrology  system  presented  unique  challenges  to  the  researchers  which  required 
significant  development. 


Accomplishments 

The  ultra-high  speed  IMACON200  cameras  have  unique  optics  that  presents  unexpected  challenges  for  their  use  with  digital 
image  correlation.  Although  not  perceivable  when  images  are  viewed  on  a  monitor,  imaging  warping  and  shearing  have  been 
measured  when  using  digital  image  correlation  with  a  single  IMACON200  camera  to  measure  unifonn,  in-plane  translations.  Through 
a  detailed  study  of  the  imaging  process  and  the  associated  errors,  we  have  developed  a  basic  understanding  for  the  source  of  the  errors. 
Working  with  a  local  company  specializing  in  optical  methods,  experimental  techniques  and  associated  software  have  been  developed 
to  characterize  the  cameras  and  improve  the  image  correlation  measurements. 

The  techniques  developed  for  the  IMACON  200  cameras  and  their  unique  optics  have  been  shown  to  work  in  all  imaging 
system  that  has  distortion.  The  scanning  electron  microscope  is  a  device  which  has  unique  distortions  which  must  be  overcome  if  it  is 
to  be  used  for  metrology  purposes.  The  techniques  developed  have  the  ability  to  allow  the  SEM  to  be  used  for  metrology  in  regards  of 
deformation  on  small  scales.  The  personal  at  ARL’s  facility  in  Aberdeen  Maryland  have  expressed  interest  in  the  use  of  this  method. 
The  paper  which  discusses  this  is  attached  to  the  final  report  since  it  has  been  accepted  but  not  yet  published. 

Manufacturing  techniques  have  been  refined  to  properly  produce  the  functional  graded  ceramic-metallic  materials.  Initial 
samples  were  produced  and  tests  of  the  samples  for  the  required  consistency  and  integrity  is  ongoing,  with  preliminary  results 
indicating  that  samples  can  be  produced  reliably.  For  example,  two  material  interface  samples  were  produced  and  delivered  to  USC 
for  initial  experimental  setup  and  verification.  Papers  have  been  published  on  the  functionally  graded  material  and  processing 
requirements  and  are  listed  in  the  paper  published  section. 

Initial  imaging  using  the  ultra-high  speed  cameras  of  the  2D  impact  of  the  two  material  interface  samples  were  obtained  using 
the  cameras  in  a  2-D  configuration.  The  experimental  parameters  developed  for  this  application  include  proper  triggering,  camera 
control  during  impact  process  and  specimen  lighting.  Due  to  large  distortions  in  the  experimental  system  the  results  were  inclusive  and 
work  will  be  ongoing  to  evaluate  these  specimens. 
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Technology  Transfer 

A  tour  of  our  laboratory-scale  impact  facility  was  given  to  Bruce  Fink,  Chief,  Materials  Division,  Army  Research  Laboratory. 
During  the  tour,  a  high  speed  impact  event  was  captured  by  a  single  IMACON200  camera,  and  the  images  converted  into  deformation 
fields  with  our  new  software.  He  expressed  his  strong  support  for  the  development  process  and  stated  that  he  “sees  tremendous 
potential  in  the  application  of  this  technology  to  a  wide  variety  of  efforts  involving  rate-dependent  properties  of  materials  and  high- 
rate  testing  of  military  structures.” 

Conversations  and  a  meeting  at  a  professional  society  conference  have  been  had  with  Dr,  Tusit  Weerasooriya,  of  the  Impact 
Physics  Branch  of  Army  Research  Laboratory.  He  has  stated  “Once  the  system  is  developed,  we  intend  to  use  it,  with  USC 
collaboration,  in  our  high  rate  lab”.  He  has  given  a  few  examples  where  they  intend  to  use  the  system  which  included 

•  Measuring  the  back  surface  deformation  of  personnel  armor  during  impact  to  validate  computational  models. 

•  To  develop  valid  test  techniques  to  measure  high  rate  behavior  of  ceramics  and  soft  polymers.  These  materials  are  extremely 
difficult  to  test  under  high  rates  using  traditional  experimental  techniques. 

•  Measure  fracture  behavior  of  penetrator  rods. 

The  experimental  work  involves  the  use  of  3-D  image  correlation  developed  by  a  small,  start-up,  high  tech  measurements 
company,  Correlated  Solutions,  Inc.  (CSI).  As  the  first  graduate  from  the  USC  Incubator  program,  the  company  is  now  located  a  few 
miles  from  the  campus  of  the  University  of  South  Carolina  and  has  been  supportive  of  the  on-going  research  activities  in  the 
experimental  mechanics  group  at  USC.  Currently,  we  are  working  closely  with  CSI  personnel  in  the  development  of  unique  software 
needed  to  control  the  ultra-high  speed  cameras  and  analyze  the  resulting  images.  Correlated  Solutions  plans  to  market  systems  using 
this  unique  software. 

We  worked  with  Dr.  Dave  McGowen  of  NASA-Langley’s  Structures  and  Materials  Division  to  make  dynamic  structural 
measurements  using  the  ultra-high  speed  3D  measurement  system  that  is  being  used  in  this  work. 
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Dr.  McNeill’s  sabbatical  at  ARL’s  facility  at  Aberdeen  Maryland 

Many  of  the  experimental  techniques  developed  for  this  work  was  of  great  interest  to  researchers  at  ARL’s  facility  at  Aberdeen 
Maryland.  Dr.  McNeill  spent  all  of  2005  at  ARL  transferring  the  technology  and  helping  with  problems  of  interest  to  the  ARL 
researchers.  Below  is  a  list  of  some  of  the  projects  and  their  results  that  was  worked  on  by  Dr.  McNeill  during  the  sabbatical. 

Project  #  1:  Determination  of  stress  strain  properties  for  green  bullet  system 

In  an  effort  to  reduce  the  potential  long-term  environmental  effects  that  may  accrue  due  to  the  use  of  lead  in  ballistic 
applications,  the  US  Army  has  begun  development  of  new  ballistic  materials  that  have  similar  material  properties  to  the  lead  currently 
being  used.  Compression  testing  is  being  performed  on  various  materials  that  have  the  same  size  and  shape  of  the  lead  core  they  are 
replacing.  The  specimens  being  tested  are  4.6  mm  in  diameter  and  6.4  mm  in  length. 

Significant  Contributions 

All  experimental  data  was  obtained  as  the  specimen  incurred  elastic  deformation.  Full-field  displacement  and  strain 
measurements  were  obtained  using  3D  Digital  Image  Correlation,  pioneered  at  the  University  of  South  Carolina.  Prior  to  the  use  of 
Digital  Image  Correlation,  the  specimens  were  instrumented  with  strain  gauges  and  tested  in  a  standard  compression  load  frame.  The 
time  required  to  instrument  each  specimen  with  strain  gauges  was  over  one  hour  due  to  the  small  size  of  the  specimens  being  studied. 
As  a  result,  only  a  few  specimens  could  be  instrumented  and  tested  per  day.  With  digital  image  correlation,  many  specimens  could  be 
speckled  in  a  matter  of  a  few  minutes  and  then  tested.  Using  3D  digital  image  correlation  it  was  possible  to  speckle  and  test 
specimens  at  a  rate  of  three  per  hour,  improving  the  through-put  of  the  experimental  program  while  also  developing  full-field  data  to 
determine  the  overall  response  of  the  green  bullet  material  under  compressive  loading. 

Project  #  2:  Small  Disk  Impact 

A  commonly  used  impact  experiment  employs  a  small  flat  round  specimen.  The  specimen  is  placed  perpendicular  to  the 
muzzle  of  a  gas  cannon  and  then  impacted  by  a  spherical  projectile.  Analysis  of  the  specimen  consists  of  measuring  the  height  of  the 
deformation  after  impact.  Since  no  strain  information  has  been  available  in  the  past,  only  qualitative  data  has  been  obtained 
previously. 

Significant  Contribution 

All  experiments  were  performed  using  3D  digital  image  correlation,  a  technique  pioneered  at  the  University  of  South 
Carolina.  Both  surfaces  of  the  specimen  were  speckled  and  images  obtained  prior  to  being  placed  in  the  gas  gun  fixture.  After 
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impact,  the  specimen  was  removed  and  images  obtained.  Then,  the  images  were  analyzed  using  digital  image  correlation  software,  so 
that  the  following  data  was  obtained 

(a)  the  complete  3D  shape  of  the  final  surface 

(b)  the  complete  surface  deformation  fields  on  both  surfaces  is  obtained  after  impact 

Figure  1  shows  the  measured  shape  of  the  surface  opposite  to  the  impact.  Figure  2  shows  the  measured  shape  of  the  impacted 
surface;  note  that  there  was  data  obtained  at  the  impact  site,  allowing  investigators  to  acquire  new  measurement  data  heretofore 
unobtainable.  Figure  3  shows  the  radial  strain  field  after  impact  on  the  surface  opposite  to  the  impact. 


Figure  1:  Radial  strain  calculated  using  strains  from  3D  image  correlation 
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Project  #  3: 


Figure  2:  Profile  of  impact  side  after  impact 


Dynamic  Deformation  Measurements  in  Kevlar  Fabric  During  Impact 
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There  has  been  a  paucity  of  data  for  validation  of  existing  computational  models  to  predict  the  deformation  process  in  Kevlar 
fabric  undergoing  projectile  impact  loading. 

Significant  Contribution 

To  obtain  the  first  set  of  dynamic  data,  two  high  speed  cameras  were  placed  in  a  stereo-vision  arrangement  and  images  were 
acquired  during  the  impact  event.  3D  digital  image  correlation,  pioneered  at  the  University  of  South  Carolina,  was  used  to  analyze  a 
series  of  images  acquired  during  impact. 

Two  sets  of  experiments  were  performed.  One  set  of  experiments  was  performed  at  ARL  using  their  gas  cannon.  A  second  set 
of  experiments  were  performed  on  an  ATC  test  range  with  projectile  velocities  exceeding  1,000  m/s. 

Typical  results  are  shown  in  Figure  4.  The  surface  shape  during  impact  is  clearly  shown  in  Fig.  4.  In  addition,  the  lagrangian 
strain  field  is  also  shown  in  color  contours  directly  on  the  deformed  surface.  As  can  be  seen  in  Fig.  4,  data  is  obtained  throughout  the 
visible  region  on  the  specimen  surface. 

Due  the  quality  of  the  data  obtained,  discussions  are  being  held  for  performing  additional  3D  deformation  measurements 
during  ballistic  impact  of  helmets  currently  being  worn  by  US  soldiers. 
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Figure  4:  One  component  of  the  normal  surface  strain  measured  in  a 

Kevlar  sheet  during  impact 
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Project  #  4:  High  speed  deformation  measurements  of  EM  gun  barrel  during 
charging  and  discharging 

The  deformation  measurement  of  the  barrel  of  the  electro-magnetic  gun  during  a  full  charging  and  discharging  cycle  would  be 
helpful  to  understand  how  the  gun  is  performing  during  operation.  However,  methods  were  not  available  that  were  capable  of 
acquiring  these  measurements. 

Significance  of  the  Contribution 

To  obtain  the  first  set  of  dynamic  data,  two  high  speed  cameras  were  placed  in  a  stereo-vision  arrangement  and  images 
acquired  of  the  EM  gun  barrel  during  a  complete  charging  and  discharging  cycle.  After  the  images  were  acquired,  3D  digital  image 
correlation,  pioneered  at  the  University  of  South  Carolina,  was  used  to  analyze  a  series  of  the  images.  Deformations  were  acquired 
as  the  EM  gun  was  subjected  to  operating  currents  ranging  from  500  kA  to  800  kA. 

Figure  5  shows  the  measured  radial  strain  on  the  barrel  of  the  EM  gun  at  the  end  of  the  cycle.  As  shown  in  Fig  5,  the  radial 
strain  field  is  not  symmetric  and  there  is  an  unexpected  large  strain  area  near  the  edge  of  the  rail;  the  location  of  the  embedded  rail 
is  between  the  two  while  lines  in  Fig  5.  The  measured  strain  fields  have  resulted  in  investigations  into  the  fabrication  techniques  used 
to  build  the  gun 


EM  Gun 
Barrel 
Axis 


Figure  5: 


Full  field  radial  strain  during  discharging 
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The  strain  for  one  point  on  the  barrel  as  a  function  of  time  for  various  charge  voltages  is  shown  in  the  Figure  6.  This  data 
suggests  that  the  EM  gun  barrel  underwent  in-elastic  deformation  during  the  800  kA  cycle. 
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Figure  6:  Hoop  strain  vs  time  for  one  point  of  EM  gun 
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Abstract 

A  novel  approach  for  correcting  both  spatial  and  drift  distortions  that  are  present  in  scanning  electron 
microscope  (SEM)  images  is  described.  Spatial  distortion  removal  is  performed  using  a  novel  methodology 
employing  a  series  of  in-plane  rigid  body  motions  and  a  generated  warping  function.  Drift  distortion  removal  is 
performed  using  multiple,  time-spaced  images  to  extract  the  time-varying  relative  displacement  field  throughout  the 
experiment. 

Results  from  numerical  simulations  clearly  demonstrate  that  the  correction  procedures  successfully  remove 
both  spatial  and  drift  distortions.  Specifically,  in  the  absence  of  intensity  noise  the  distortion  removal  methods 
consistently  give  excellent  results  with  errors  on  the  order  of  ±0.01  pixels. 

Results  from  the  rigid  body  motion  and  tensile  loading  experiments  at  200X  indicate  that,  after  correction  for 
distortions,  (a)  the  displacements  have  nearly  random  variability  with  standard  deviation  of  0.02  pixels;  (b)  the 
measured  strain  fields  are  unbiased  and  in  excellent  agreement  with  previous  full-field  experimental  data  obtained 
with  optical  illumination;  (c)  the  strain  field  variability  is  on  the  order  of  60  microstrain  in  all  components  with  a  spatial 
resolution  on  the  order  of  25  pixels. 

Taken  together,  the  analytical,  computational  and  experimental  studies  clearly  show  that  the  correction 
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procedures  successfully  remove  both  spatial  and  drift  distortions  while  retaining  excellent  spatial  resolution, 
confirming  that  the  SEM-based  method  can  be  used  for  both  micro-material  and  nano-material  characterization  in 
either  the  elastic  or  elastic-plastic  deformation  regimes. 
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1.  Introduction 

Computer  vision  and  digital  image  correlation  (DIC)  have  been  applied  to  the  study  of  in-plane  material  behavior  at 
the  macro-scale  [1'4].  For  nominally  planar  objects  undergoing  in-plane  motions,  the  2D-DIC  technique  became  one  of  the 
preferred  optical  methods  for  measurements  in  experimental  mechanics.  For  general  surface  shapes  with  arbitrary 
motions,  3D-DIC  (stereo-vision)  employing  at  least  two  views  to  recording  object  positions  has  been  shown  to  have 
general  capability  to  measure  3D  shape  and  full-field  3D  motions  [5'13].  Nowadays,  2D  and  3D  methods  using  computer 
vision  are  widely  used  in  many  applications,  and  commercial  software  is  available1141. 

A  first  step  to  access  a  reduced-length  scale  of  measurement  is  to  apply  the  DIC  method  with  images  obtained  using 
high-magnification  optical  systems  for  2D  or  3D  measurements  [15'17].  pew  authors  have  investigated  the  problem  of  the 
accurate  calibration  of  micro-scale  imaging  systems,  including  the  determination  and  correction  of  the  underlying 
distortions  in  the  measurement  process.  One  reason  may  be  the  obvious  complexity  of  high-magnification  imaging 
systems  that  weakens  the  underlying  assumptions  commonly  used  in  parametric  distortion  models  that  correct  simple 
lens  systems  such  as  digital  cameras [18, 19]. 

Recently,  Schreier  et  al.  proposed  a  new  methodology  to  calibrate  accurately  any  imaging  sensor  by  correcting  a 
priori  for  the  distortion  using  a  non-parametric  model [13].  The  a  priori  correction  of  distortion  transforms  the  imaging  sensor 
into  a  virtual,  distortion-free  sensor  plane  using  unknown  arbitrary  rigid  body  motions  of  a  gridless  planar  target.  As 
opposed  to  classical  calibration  techniques  relying  on  a  dedicated  target  marked  with  fiducial  points,  this  approach  can  be 
applied  using  any  randomly  textured  planar  object  (i.e.,  a  “speckle  pattern”).  This  type  of  distortion  depends  only  on  the 
pixel  location  in  the  image  and  is  designated  as  spatially-varying  distortion  or  spatial  distortion  in  this  work. 

Due  to  the  nature  of  light,  optical  imaging  systems  are  limited  to  a  maximum  resolution  that  corresponds  to  a 
magnification  of  «  1,000X.  For  higher  spatial  resolution  (and  also  smaller  displacements),  imaging  systems  based  on 
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electron  microscopy  (such  as  SEM  and  TEM)  are  employed.  Since  the  physics  of  electron  microscopy  [20,  211  is  quite 
different  from  optical  microscopy,  a  new  model  and  calibration  process  is  necessary  compared  to  the  classical  approach. 
For  example,  SEM  systems  have  both  spatial  distortion  and  also  a  time-varying  distortion  or  drift  distortion.  Drift  distortion 
occurs  in  an  SEM  imaging  system  since  the  electron  beam  scanning  process  oftentimes  results  in  a  non-uniform 
“apparent”  displacement  field  across  the  specimen;  relative  motions  occur  between  pixel  locations  that  have  been  shown 
to  change  with  time.  The  cumulative  effect  of  both  spatial  and  drift  distortions  can  introduce  large  image  displacements 
and  substantial  image  deformations. 

To  correct  for  the  combined  effect  of  spatial  and  drift  distortion,  different  projection  models  have  been  established 
but  no  work  has  provided  both  a  satisfactory  model  and  a  correction  method  for  such  image  distortions.  In  fact,  most 
papers  and  even  commercial  SEM  measurement  systems  simply  ignore  these  effects  and  consider  a  pure  projection 
model t22'28].  A  few  authors  taking  into  account  distortion  consider  only  a  parametric  model [29'311  and  neglect  the  effect  of 
drift  [32'37]. 

In  this  paper,  the  procedure  we  have  developed  for  quantifying  time-varying  distortion  (Section  2.1)  and  spatially- 
varying  distortion  (Section  2.2)  is  presented,  along  with  an  overall  procedure  for  implementing  the  distortion  determination 
procedure  (Section  2.3).  Section  3  describes  computer  simulations  performed  to  demonstrate  the  potential  accuracy  of 
the  method  when  employed  to  identify  and  remove  typical  spatial  and  drift  distortion  fields.  Section  4  describes  the 
translation  and  uniaxial  strain  experiments  performed  in  an  SEM  to  demonstrate  the  accuracy  of  the  methods  developed. 

2.  Image  Distortions  in  an  SEM 

It  is  assumed  that  there  are  two  independent  distortion  functions.  First,  Ddr(t)  is  defined  as  the  drift  distortion  function, 
where  t  is  the  scan  time  and  t  e  [0,°°).  Second,  Dsp(x,y)  is  defined  as  the  spatial  distortion  function,  where  (x,y)  is  a  pixel 
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location  on  the  image  plane. 

A  position  is  distorted  due  to  both  drift  and  spatial  distortion  in  an  SEM  image.  Since  distorted  positions  are  the 
measurable  quantities  in  experiments,  the  inverse  functions  for  both  drift  and  spatial  distortion  are  the  quantities  that  are 
determined  during  the  calibration  process4. 

2.1.  Time-Varying  Distortion 

Time-varying  distortion  or  drift  distortion  occurs  throughout  the  scan  process,  introducing  artificial  relative 
displacements  as  the  image  is  formed,  resulting  in  errors  when  computing  strains  using  the  drift-distorted  images;  we 
have  observed  relative  displacements  within  a  single  image  of  up  to  6  pixels  at  10,000X  though  relative  drift  of  0.80  pixels 
across  an  image  is  more  common.  Because  of  its  non-stationary  nature  (magnitude  and  direction  are  a  function  of  time), 
this  effect  cannot  be  removed  when  performing  spatial  distortion  correction. 

To  quantify  the  drift  distortion  at  each  pixel  throughout  an  experiment,  a  novel  drift  correction  procedure  has  been 
developed  and  is  described  in  the  following  sections. 

2. 1.1. SEM  Scanning  Basics 

A  SEM  image  is  generated  through  a  raster-scanning  process.  Each  pixel  requires  a  dwell  time,  to  ,  to  define  the 
“intensity”  of  image  at  that  location,  so  that  the  time  required  to  scan  an  entire  row,  tR,  and  an  entire  frame,  tF,  of  image 
are  given  by 

4  Though  the  focus  of  this  study  is  2D-DIC  with  a  single  view,  the  procedure  should  be  applied  to  each  camera  (view)  in  a  multi-camera  system 
to  remove  distortions  from  individual  views.  Thus,  the  procedure  is  equally  important  when  3D  reconstruction  are  performed  using  SEM 
images. 


15/44 


tR  -  Wtp  +  tj 

tF  =  HtR  =  HWtD  +  Htj 


(1) 

(2) 

where  W  is  the  number  of  columns  in  the  image,  H  is  the  number  of  rows  in  the  image  and  tj  is  the  time  delay  to 
reposition  and  stabilize  the  electron  beam  prior  to  initiating  the  next  row  scan.  Since  the  (x,y)  position  in  the  image  is  in 
direct  correspondence  with  the  scan  time,  one  can  write 

t(x,y)  =  xtD  +  ytR  (3) 

where  0  <  x  <  W-1  and  0  <  y  <  H-1 . 

Eq.  (3)  implies  that  two  pixel  locations,  and  hence  two  distinct  scan  times,  will  generally  experience  much  different 
drift  functions.  Such  differences  have  been  observed  for  consecutive  rows  in  an  image,  with  a  clearly-defined  shift  in  drift 
measured  when  moving  from  the  last  pixel  in  one  row  to  the  first  pixel  in  the  next  row. 

In  order  to  correct  SEM  images  for  drift,  three  approaches  have  been  investigated:  (a)  global  model  of  the  drift 
function  developed  with  the  data  from  all  images  to  correct  for  the  drift  at  (x,y)  positions  in  every  image;  (b)  local  time- 
based  models,  Ddr(t),  that  are  different  for  each  image  and  (c)  local  spatial-based  models,  Ddr(x,y),  that  are  also  different 
for  each  image.  It  is  noted  that  (b)  and  (c)  imply  that  the  drift  distortion  will  only  be  known  locally  and  must  be  re-estimated 
when  a  new  image  is  acquired.  This  “incremental”  procedure  to  determine  drift  distortion  correction  requires  that 
corrections  be  made  throughout  the  experimental  process,  and  is  in  contrast  with  spatial  distortion,  which  is  determined 
during  the  calibration  phase  and  does  not  require  re-estimation  during  the  post-calibration  phase  of  an  experiment. 

Preliminary  experiments  and  comparison  of  experimental  observations  to  model  predictions  clearly  indicate  that  the 
global  model  cannot  adequately  represent  experimental  observations  whereas  both  of  the  local  model  predictions  are  in 
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good  agreement  with  experimental  measurements.  In  the  following  sections,  only  the  two  local  models  will  be  discussed 
in  more  detail. 

2.1.2.Time-based  Local  Model 

The  drift  distortion  function  is  expressed  as  follows 

Ddr(t)  =  [  5x(t),  5y(t)  ],  units  in  pixels  (4) 

where  t  is  given  in  Eq.  (3)  for  a  point  of  interest  in  the  image,  (x,y). 

Consider  a  object  point  P  located  at  position  (x,y)  in  the  image  plane  at  time  t.  At  time  t+A,  the  image  of  the  object 
point  is  located  at  a  new  position,  P’,  due  to  drift. 

The  difference  in  positions  for  the  same  object  point  is  the  disparity,  dispdr,  and  can  be  written  in  the  form 
dispdr  =  Ddr(t+A)  -  Ddr(t)  (5) 

with 

t  =  xtD  +  ytR 
A  —  tp  +  Atn  +  tdisp 

tdisp  =  [  5x(t+A)  -  5x(t)  ]tD  +  [  5y(t+A)  -  5y(t)  ]tR 

where  Atn  is  the  delay  time  between  one  image  scan  and  the  next  image  scan.  Fig.  1  shows  a  schematic  of  the  process 
mathematically  described  by  Eqs.  (4-5),  with  the  experimentally-observed  shift  between  images  clearly  visible  since  the 
drift  for  position  (W-1 , H— 1 )  is  not  similar  to  the  drift  for  position  (0,0)  in  the  next  image. 

From  the  drift  evolution  shown  in  Fig.  1,  different  functional  forms  can  be  assumed  to  represent  the  local  drift.  Fig.  1 
shows  an  estimated  linear  fit  to  the  drift  data  for  three  images.  The  linear  fit  will  smooth  regions  with  high  gradients  in  drift, 
especially  those  positions  where  step  changes  may  occur  during  the  image  acquisition  process. 
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To  determine  the  fitting  parameters  for  the  assumed  form  of  the  drift  distortion  function,  two  consecutive  images  are 
taken  at  each  step  in  the  experiment.  Assuming  that  only  drift  occurs  between  the  two  images  in  each  consecutive  image 
pair,  the  drift  velocity,  Vdr(t),  at  each  position  is  determined  by  the  disparity  in  each  pair  divided  by  the  time  increments. 

vdr(t+A/2)  =  { [  6x(t+A)  -  5x(t)  ]/A,  [  5y(t+A)  -  5y(t)  ]/A }  (6) 

where  A  is  defined  for  pixel  (x,y)  in  Eq.  (5)  and  the  disparity  values  are  obtained  in  the  sensor  plane. 

Advantages  of  the  time-based  approach  are  its  direct  relationship  to  the  actual  scanning  process  and  the  relative 
simplicity  of  the  data-fitting  process  for  time-streamed  disparity  data.  Disadvantages  are  as  follows.  First,  there  are  gaps 
in  the  time-history  due  both  to  the  re-positioning  of  the  electron  beam  after  each  row  scan  and  to  the  time  delay  between 
acquisitions  of  images.  Second,  the  drift  data  in  actual  experiments  exhibits  step  changes  in  magnitude  during  the  time- 
gaps,  resulting  in  difficulties  in  extracting  accurate  results  near  these  positions.  Third,  finite  sized  subsets  in  image 
matching  process  will  introduce  averaging  effects  across  the  scanned  rows,  introducing  errors  in  the  time-based 
approach. 

2. 1.3. Hybrid  Spatial-Temporal  (HST)  Local  Model 

As  shown  in  Eqs.  (3-5),  the  drift  function  can  be  written  in  terms  of  spatial  variables,  the  integer  pixel  position  (X|,Yj). 
The  elapsed  time  during  the  whole  imaging  process  for  calibration  can  be  written  as 
t  =  X,  tD  +  Yj  tR  +  (n— 1  )(tF  +  Atn),  n=1 ,2,...N  (7) 

where  N  is  the  total  number  of  images  for  the  drift  function  estimation. 

In  this  form,  one  can  write  Ddr(X|,Yj)  as  the  total  drift  function  and  dispdr(Xi,Yj)  as  the  disparity.  The  overall  imaging 
process  for  the  HST  model  is: 

o  A  specimen  is  patterned  so  that  the  SEM  images  are  appropriate  for  image  correlation. 
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o  A  sequence  of  image  pairs  is  acquired.  They  are  numbered  as  (1,2),  (3,4),...,(n-1,n)...,(N-1,N),  where  drift  occurs 
between  image  n-1  and  n  in  an  image  pair. 

o  Between  the  acquisition  of  each  image  pair,  rigid  body  motions  are  applied  (these  will  be  used  to  quantify  the 
spatial  distortion  field). 

o  For  each  image  pair,  2D-DIC  is  performed5  to  determine  the  drift  disparity,  dispdr(X|,Yj),  for  each  fixed  position 
(X|,Yj)  in  the  sensor  plane. 

o  A  B-spline  function  is  used  to  obtain  a  least-square  best  fit  surface  to  the  components  of  dispdr(X|,Yj)  for  each 
image  pair. 

o  With  known  values  for  to  ,  tR  and  tF,  along  with  recorded  values  for  Atn  between  each  image,  Eq.  (7)  can  be  used 
to  determine  the  elapsed  time  associated  with  position  (X|,Yj)  in  each  image.  Using  Eq.  (6),  the  velocity  of  drift  for 
each  position,  Vdr(X|,Yj;t+A/2),  can  be  estimated  using  a  central  finite  difference  form. 

o  At  each  position  (X|,Yj),  Vdr(X|,Yj;t)  is  fitted  with  a  B-spline  function  in  time. 

o  Using  the  initial  condition  Ddr(0,0)=(0,0)  at  t=0,  the  B-spline  approximation  for  vdr(X|,Yj;t)  can  be  integrated  over 
time  to  determine  the  drift  for  each  position  (X|,Yj)  in  the  nth  image6,  Ddr,n(Xi,Yj),  n=1 ,3,...,N-1 . 

o  Ddr,n(Xi,Yj)  is  used  to  correct  images  n=1 ,3,..., N-1 ,  for  the  drift  distortion. 

Mathematically,  the  determination  of  the  time  dependent  drift  function  for  each  fixed  position  (X|,  Yj)  is  obtained  via 

5  When  using  spatially-based  drift  correction,  the  subset  size  should  be  as  small  as  possible  without  introducing  non-uniqueness  in  the  matching 
process,  thereby  minimizing  local  averaging. 

6  For  the  initial  image  pair,  variations  in  disparity  across  the  image  requires  a  piecewise  integration  over  the  disparity  field  between  (0,0)  and 
(X|,Yj)  so  that  the  drift  at  any  position  (X|,Yj)  can  be  determined.  A  less  accurate  approach  would  be  to  integrate  the  estimated  dispdr(X|,Yj;t) 
from  t=0  to  the  time  corresponding  to  (X|,Yj). 
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(8) 


minimization  of  the  following  functional 

«>  =  E|va,(X1,YJ;tn)-va“(Xl,YJ;a|,tr,)|2 

n=1 

where  aj  are  the  parameters  to  be  fitted  and  tn  is  the  elapsed  time  associated  with  each  drift  disparity  at  the  integer  pixel 
position  (X|,Yj).  To  obtain  the  total  drift  function  at  each  position 

Dar(X„YJ;t)  =  }[vdr"(Xl,YJ,a1,t)]dt,  where  Ddr(0,0;0)  =  0  (9) 

0 

Once  the  drift  distortion  is  known  as  a  function  of  time  for  each  pixel  position  in  each  image,  this  field  is  fitted  by  a 
piecewise  B-spline  to  provide  a  functional  form  for  each  image,  Ddrn(x,y),  where  (x,y)  is  any  real-valued  position  in  the 
image. 

The  HST  model  described  above  offers  one  major  advantage:  the  spatial  description  of  drift  within  each  image 
minimizes  the  effects  of  time  gaps  on  local  estimates  of  the  drift  function. 

2.2.  Spatially-varying  Distortions 

In  simple  lens  systems,  such  as  used  in  a  digital  camera,  spatially-varying  distortion  or  spatial  distortion  (a.k.a., 
image  distortion)  is  a  well-known  problem.  The  commonly  used  method  for  modeling  such  imaging  systems  assumes  that 
the  distortions  are  not  a  function  of  time.  Classical  models  used  to  estimate  the  spatial  distortions  are  parametric  in 
nature  [38];  typical  model  forms  include  radial  distortion,  de-centering  distortion,  prismatic  distortion  and  tangential 
distortion. 

Since  the  SEM  imaging  process  is  based  upon  the  interaction  between  atoms  of  the  observed  specimen  and  an 
electron  beam,  as  well  as  scanning  and  focusing  processes  that  employ  electro-magnetic  principles  to  perform  the 


20/44 


required  functions,  it  is  clear  that  the  pre-specified  classical  distortions  are  likely  to  be  ineffective  in  quantifying  arbitrary 
aberrations  or  unknown  (but  deterministic)  distortions  in  a  complex  imaging  system. 

To  address  this  issue,  non-parametric  forms  for  distortion  models  are  preferred.  Two  of  the  earliest  such  models  are 
proposed  by  Peuchot [39]  and  Brand  et  al. [401,  where  cross-targets  with  known  spatial  separation  are  used  for  the  distortion 
estimation.  Since  such  targets  are  likely  to  be  difficult  to  realize  at  the  micro-  or  nano-scale  for  the  SEM,  the  work  of 
Schreier  et  al.  using  a  non-parametric  distortion  correction  approach  and  a  speckle-patterned  calibration  target [131  is  the 
most  readily  adaptable  for  distortion  correction  in  an  SEM,  and  is  employed  in  previous  work[41]  as  well  as  in  this  study. 

As  outlined  by  Schreier  et  al [13],  the  process  uses  images  acquired  during  translation  along  two  orthogonal  directions 
and  B-splines  or  other  general  forms  to  determine  full-field,  spatial  distortions.  Here,  it  is  assumed  that  the  procedures 
outlined  in  Section  2.1 .3  have  been  employed  to  remove  drift  distortions  from  all  of  the  translated  images. 

2.3.  Overall  Process  for  Distortions  Removal 

The  procedure  used  in  this  study  to  extract  drift  and  spatial  distortions  is  shown  as  a  flow  chart  in  Fig.  2.  The  pairs  of 
images  are  separately  correlated  throughout  the  calibration  phase  and  measurement  phase  to  obtain  the  “drift”  at  selected 
pixel  positions.  The  resulting  drift  disparity  fields  are  used  to  obtain  the  drift  velocity  at  selected  pixel  positions  throughout 
the  field  of  view.  Using  the  velocity  field  as  estimates  for  the  local  time  derivatives,  B-spline  fits  are  obtained  for  the  drift 
distortion. 

After  removing  the  estimated  drift  distortion  from  each  pixel  position,  the  spatial  distortion  field  and  rigid  body 
motions  during  the  calibration  phase  are  determined  through  an  optimization  procedure.  Relaxation  principles  are 
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employed  during  the  iteration  process  to  verify  that  the  estimated  drift  and  spatial  distortion  fields  are  converged7. 

3.  Numerical  Simulations 

Simulations  are  performed  for  both  calibration  phase  and  measurement  phase8  in  consistent  with  actual 
experimental  conditions.  The  image  correlation  process  is  not  simulated.  Both  drift  and  spatial  distortion  functions  are 
assumed  to  have  pre-defined  forms.  Using  these  functional  forms,  the  distorted  positions  of  a  finite  number  of  locations, 
(x,y)k,  are  determined  for  each  image  n.  The  process  is  repeated  at  N  times,  so  that  the  distorted  positions  (x,y)kn, 
k=1,2,...,K,  n=1,2,...,N,  are  determined9.  This  data  forms  the  basis  for  the  simulation  process,  with  the  long  term  goal  of 
demonstrating  the  accuracy  and  robustness  of  the  approach  for  SEM  image  analysis. 

3.1 .  Simulation  for  the  Calibration 

With  At=120s,  tD=  1 0'4s,  tR=1 ,071x10'2s  and  an  imaging  array  size  of  1024x884,  the  positions  of  image  points  at  30 
specific  times  (corresponding  to  15  image  pairs  in  an  experiment)  are  generated  over  a  total  time  of  128  minutes. 
Between  each  image  pair,  the  effect  of  a  cross-shaped  translation  is  included  in  the  position  of  each  point.  To  distort  the 
position  of  each  image  point,  the  drift  distortion  function  Ddr(t)  is  assumed  to  have  a  quadratic  form.  The  spatial  distortion 
function  is  assumed  to  have  the  form  of  a  combination  of  cosine  wave  and  quadratic  surface.  Fig.  3  shows  the  drift  and 

7  SEM  images  and  the  corresponding  disparity  maps  have  considerable  electronic  and  measurement  noise.  Thus,  the  HST  local  drift  correction 
model  should  employ  a  reduced  order  for  the  B-spline  fit  to  the  disparity  data  so  that  smoothing  of  the  data  is  performed  during  the  fitting 
process  so  that  oscillations  in  the  estimated  velocity  field  are  minimized. 

8  The  procedure  whereby  the  drift  distortion  is  computed  separately  for  the  calibration  and  measurement  phases  is  used  in  practice  to  minimize 

the  effects  of  specimen  shifts  during  the  initial  loading  process.  However,  in  principle  the  process  can  be  continuous. 
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spatial  distortion  functions.  Thus,  with  the  inclusion  of  random  error,  the  distorted  positions  of  an  image  point  (x,y)  in 
image  n  are  written 

P’(x,y;t)  =  P(x,y;t)  +  Ddr(x(t),y(t))  +  Dsp(x,y)  +  G(x(t),y(t))  (10) 

where  G(x,y)  is  a  Gaussian  error  function  with  mean  value  0  and  a  pre-defined  standard  deviation,  varying  with  time  and 
by  Eq.  (3)  spatial  position  (x,y). 

Instead  of  performing  correlation  to  extract  the  disparity  maps,  a  total  of  150x120=18,000  distorted  positions  with  a 
spacing  of  5  pixels  and  a  “beginning”  image  position  located  at  (101,101)  are  calculated  using  Eq.  (10)  with  known  vectors 
on  the  right  hand  side. 

To  obtain  the  drift  component  numerically,  all  15  drift  disparity  maps  are  computed  at  each  pixel  location  and  Eqs. 
(6-10)  are  applied  to  extract  Ddr(x,y).  Fig.  4  shows  a  direct  comparison  between  the  computed  and  the  input  drift  at 
position  (51 1,441).  Fig.  5  shows  the  spatial  distribution  for  the  difference  between  the  computed  and  the  input  vertical  drift 
for  the  image  acquired  at  t=36  minute. 

After  correction  to  remove  the  effects  of  drift,  the  disparity  maps  are  obtained  by  subtracting  the  positions  of  points  in 
the  reference  image,  image  1,  from  the  positions  of  matching  points  in  all  other  odd-numbered  images.  The  resulting 
disparity  maps  have  the  form 

dispsp(x,y)  =  Dsp(x+u,y+v)  -  Dsp(x,y)  +  G(x,y)  (1 1 ) 

where  (u,v)  is  the  orthogonal  rigid  body  motion  of  other  images  relative  to  image  1 . 

By  incorporating  the  computed  disparity  maps  into  Eq.  (10),  the  procedure  described  in  [13]  is  employed  to  perform 


9  In  practice,  7  to  1 1  pairs  of  images  are  acquired  during  the  calibration  phase;  the  only  requirement  is  that  several  translations  in  two  orthogonal 
directions  be  performed.  The  number  of  pairs  of  images  during  the  measurement  phase  will  vary  with  the  number  of  strain  increments;  for 
better  estimation  of  B-Spline  function  it  should  be  more  than  6  pairs  of  images  acquired. 
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least-square  bundle-adjustment  optimization  and  determine  all  translations  and  all  parameters  in  the  spatial  distortion 
function.  Fig.  6  shows  the  difference  between  the  computed  and  the  input  spatial  distortion. 

As  a  final  check  on  the  accuracy  of  the  drift  and  spatial  distortion  correction  method,  the  residual  strain  fields  (with 
and  without  Gaussian  noise  in  the  displacement  values)  are  computed  for  all  of  the  images  within  the  calibration  phase. 
Without  Gaussian  noise  in  the  displacement  components,  the  computed  strains  are  less  than  1  microstrain  throughout  the 
entire  sequence.  With  Gaussian  noise  having  a  standard  deviation  *  0.025  pixels  in  each  displacement  component,  all 
strains  have  an  average  strain  between  ±60  microstrain  and  a  standard  deviation  «  50  microstrain. 

3.2.  Discussion 

The  simulations  assumed  a  total  drift  of  10-20  pixels  over  2  hours.  Thus,  between  images  in  the  sequence,  the  drift 
is  relatively  small  and  the  local  drift  velocities  can  have  considerable  oscillation  due  to  electronic  noise.  Even  so,  the 
simulations  confirmed  that  the  method  proposed  will  give  good  overall  accuracy,  even  in  the  presence  of  substantial 
Gaussian  noise  in  the  measurements,  when  combining  both  drift  and  spatial  distortion  correction. 

Since  the  drift  is  relatively  small  between  images,  it  may  appear  that  one  can  simply  ignore  this  phenomenon. 
However,  our  simulations  indicate  that  ignoring  or  incorrectly  estimating  drift  distortion  will  introduce  substantial  errors  in 
the  spatial  distortion  correction.  These  errors  will  introduce  large  residual  strains  («  1,000  microstrain)  that  cannot  be 
removed  from  the  data. 

4.  Experiment  Validation  for  Calibration  and  Strain  Measurement 
4.1.  Experimental  Setup 

Consistent  with  results  obtained  in  the  previous  sections,  all  SEM  imaging  in  this  paper  is  performed  using  an  FEI 
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Quanta-200  SEM  in  high  vacuum  mode.  For  the  electron  beam,  (a)  the  accelerating  voltage  is  30KV,  (b)  beam  spot  size  is 
3,  (c)  tD=1 0'4s,  (d)  tR=1 .071x1 0'2s.  The  CCD  sensor  array  has  a  size  of  1024x884.  These  data  corresponds  to  tF=94.67s. 
All  images  are  obtained  using  the  BSE  detector  with  a  working  distance  14.3mm.  For  validation  purposes,  a  magnification 
of  200X  is  used,  which  corresponds  to  «  1 ,25pm  per  pixel.  The  intensity  value  is  stored  in  8-bits  format  in  the  file. 

To  apply  a  micro-scale  pattern  to  the  specimen,  the  procedure  described  in  previous  publications  [42'461  is  used.  For 
comparison  to  the  vision-based  measurements,  a  single  strain  gage  is  aligned  with  the  loading  direction.  All  experiments 
are  performed  with  the  specimen  installed  in  a  tensile  loading  frame  and  mounted  inside  the  SEM  chamber. 

4.2  Experiments 

The  first  experiment  corresponds  to  the  calibration  phase  and  requires  a  series  of  translations  in  two  orthogonal 
directions  performed  manually  via  external  motion  controls.  The  second  experiment  adds  a  series  of  strained  images  that 
correspond  to  the  measurement  phase;  uniaxial  loading  is  performed  in  load  control  using  a  Labview  control  and  data 
acquisition  program. 

After  completing  each  experiment  and  acquiring  all  images,  the  disparity  maps  are  obtained  by  2D-DIC  using  a 
specially  modified  version  of  commercial  software  VIC2D  [14].  All  correlations  are  performed  at  150x120  image  locations 
using  a  43x43  subset  size  and  spacing  between  subset  centers  of  5  pixels,  the  first  subset  center  is  at  (101,101 ).  Thus,  all 
disparity  maps  contain  18,000  points. 

The  procedures  outlined  in  Fig.  2  are  employed  to  determine  the  drift  and  spatial  distortion  functions.  To  obtain  strain 
fields  from  the  corrected  displacement  data,  local  quadratic  fits  to  the  displacement  components  are  performed.10. 

10  In  this  case,  the  spatial  resolution  is  governed  by  the  subset  size.  An  estimate  for  the  spatial  resolution  is  =  25  pixels. 
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4.3.  Calibration  Phase-Drift  and  Spatial  Distortions  Removal 

A  cross-shaped  motion  path  is  performed  with  16  motions.  Before  beginning  the  motion  sequence,  two  consecutive 
images  are  acquired  as  “references”,  resulting  in  a  total  of  17  pairs  of  images  for  calibration. 

Fig.  7  shows  the  measured  drift  velocity  in  y-direction  (dv/dt)  at  location  (236,211),  along  with  the  best-fit  quadratic 
function.  Though  oscillations  are  clearly  present  (noise  is  on  order  of  ±0.02  pixels),  they  are  small  and  only  visible  due  to 
the  relatively  small  drift  velocity  of  0.01  pixels/minute  measured  in  the  SEM  at  this  magnification. 

The  measured  spatial  distortion  functions  are  presented  in  Fig.  8.  Here,  it  is  clear  that  the  spatial  distortions  are 
much  larger  than  the  expected  drift  distortions,  with  spatial  distortion  corrections  up  to  1  pixel  in  x-direction  and  2  pixels  in 
y-direction. 

The  strain  fields  are  also  computed  in  each  of  the  16  calibration  images.  Results  indicate  that 
o  the  uncorrected  exx,  eyy  and  exy  have  average  values  approaching  200  microstrain  and  a  standard  deviation  of  100 
microstrain. 

o  the  drift-corrected  exx,  eyy  and  exy  have  average  values  between  100  and  200  microstrain  and  a  standard  deviation 
of  80  microstrain. 

o  the  spatial  distortion-corrected  exx,  eyy  and  exy  have  average  values  near  50  microstrain  and  a  standard  deviation  of 
60  microstrain. 

o  the  fully-corrected  exx,  eyy  and  exy  have  average  values  near  zero  and  a  standard  deviation  of  60  microstrain. 

4.4  Measurement  Phase-Application  of  Uniaxial  Tension 
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After  completing  the  calibration  phase,  the  specimen  is  subjected  to  increasing  uniaxial  loading  and  images  are 
acquired  at  6  different  load  levels.  To  minimize  the  potential  for  unwanted  image  motion,  after  increasing  loading  on  the 
specimen  each  time,  a  hold  time  of  30s  is  maintained  prior  to  acquiring  additional  images. 

To  obtain  the  additional  drift  components  that  accumulate  during  the  strain  portion  of  the  experiment,  the  procedures 
outlined  Section  2.1.3  are  employed  to  extract  the  drift  corrections  during  the  loading  sequence.  Both  the  spatial  distortion 
and  also  the  new  local  drift  functions  are  used  to  correct  all  spatial  positions  during  the  loading  process.  The  disparity 
maps  are  acquired  across  corrected  odd-numbered  images  to  determine  the  residual  deformation  field. 

Table  1  presents  a  summary  of  the  average  data  across  the  field,  along  with  estimates  for  the  variability  of  each 
strain  component.  The  data  in  Table  1  indicates  that  the  shear  strain  is  nearly  zero  throughout  the  experiment  and  the 
standard  deviation  within  the  field  is  less  than  1 14  microstrain  for  all  components. 

Fig.  9  presents  the  average  strain  data  for  exx  (transverse  to  loading  direction)  and  eyy  (along  loading  direction)  at  all 
loading  levels,  along  with  a  direct  comparison  to  the  independently  measured  eyy  data  via  strain  gage.  As  shown  in  Fig  9, 
the  average  data  from  strain  gage  and  our  fully-corrected  data  are  in  excellent  agreement,  confirming  that  the  Young’s 
modulus  obtained  from  the  strain  gage  measurements  and  the  fully-corrected  data  are  nearly  the  same,  70GPa.  This 
value  is  consistent  with  a  wide  range  of  literature  values  for  aluminum  alloys. 

Also  shown  in  Fig.  9  are  the  strain  values  that  would  be  predicted  using  (a)  no  corrections  and  (b)  drift  correction 
only.  It  is  clear  that  fully-corrected  images  are  necessary  for  accurate  strain  estimation  so  that  reasonable  estimates  can 
be  made  for  elastic  material  response. 

Finally,  since  Poisson’s  ratio,  v,  is  defined  by  v  =  -£xx  /  £yy,  the  average  measured  value  using  fully  corrected  image 
correlation  images  is  0.33.  This  is  in  good  agreement  with  handbook  data,  which  has  a  range  0.29  <  v  <  0.34. 
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4.5.  Discussion 

Algorithms  implementing  both  spatial  and  drift  distortion  correction  method  have  been  shown  to  be  effective  for 
removing  distortions  from  SEM  images.  Thus,  the  data  clearly  indicates  that  basic  elastic  material  properties  can  be 
reasonably  quantified  using  digital  image  correlation  with  corrected  SEM  images.  Furthermore,  even  at  low  magnification, 
both  corrections  are  essential  for  accurate  measurement  of  elastic  response. 

The  methods  outlined  in  this  work  have  been  extended  by  the  authors  to  make  quantitative,  high  accuracy  3D 
displacement  and  shape  measurements  with  an 

SEM  [47],  the  complete  development  will  be  the  subject  of  forthcoming  articles.  For  such  measurements,  the  equations  of 
3D  computer  vision  replace  the  relatively  simple  2D  equations  given  in  Section  2,  and  the  calibration  phase  described 
previously  [13]  is  employed  where  the  specimen  is  rotated  about  its  eucentric  axis  to  obtain  multiple  views,  and  image 
correlation  is  employed  to  locate  a  dense  set  of  matching  positions  in  a  component.  By  several  measured  3D  profile  of 
some  objects  using  stereo  images  obtained  in  our  SEM,  though  no  quantitative  comparison  is  provided,  the  results  are 
encouraging  and  provide  confidence  that,  when  completed,  the  method  will  be  both  accurate  and  robust  in  such 
applications. 

5.  Concluding  Remarks 

The  novel  method  outlined  in  this  work  relies  on  a  combination  of  a  priori  drift  and  spatial  distortion  correction  so  that 
accurate  elastic  and  elastic-plastic  deformation  measurements  can  be  obtained  using  SEM  images;  both  corrections  are 
essential  to  obtain  accurate  deformation  measurements  throughout  the  field  of  view. 

In  sharp  contrast  with  the  approach  of  early  SEM  measurements,  where  the  investigators  simply  accepted  the 
accuracy  obtainable  and  successfully  performed  their  studies  for  important  problems  amenable  to  such  limitations,  this 
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work  presents  a  general  approach  that  successfully  extends  the  range  of  measurements  obtainable  in  an  SEM  to  the 
small  deformation  (elastic)  regime  so  that  full  elastic-plastic  deformation  studies  can  be  performed  in  an  SEM. 

Simulation  results  have  shown  that  typical  drift  processes  in  an  SEM  can  be  adequately  reconstructed  using  local 
drift  velocity  measurements.  However,  if  higher  gradients  in  drift  are  present  during  the  early  stages  of  image  acquisition, 
the  simulations  also  show  that  image  acquisition  time  should  be  reduced  and  additional  images  need  to  be  acquired 
during  this  period  for  accurate  drift  reconstruction,  a  situation  that  may  not  be  feasible  with  a  given  microscope.  In 
practice,  image  acquisition  should  be  conducted  15-30  minutes  after  the  first  image  scan  is  initiated  so  that  the  gradients 
in  drift  are  reduced  to  a  more  manageable  level. 
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Fig.  1  Schematic  representation  for  noiseless  vertical  drift  as  a  function  of 
time  for  three  images,  along  with  definition  of  disparity.  Dashed  line 
represents  a  best  linear  fit. 


Fig.  2  Overall  procedure  employed  for  correcting  distortion  in 
an  SEM.  Relaxation  methods  are  implemented  during 
process  to  improve  convergence. 
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Fig.  4  Comparison  of  the  computed  and  the  input  drift  at  pixel  (511,441)  throughout  calibration 
process  (a)  without  Gaussian  noise  and  (b)  with  Gaussian  noise  range  of  ±0.1  pixels  added 
to  the  measurements. 
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Fig.  5:  Spatial  distribution  of  difference  between 
computed  and  input  vertical  drift  for  the  image 
acquired  at  t=36  minute. 
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noise  of  ±0.1  pixels  applied  to  input  data.  Residuals  are  less  than  ±0.02  pixels 
throughout  the  field. 
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Fig.  7:  Measured  drift  velocity  as  a  function  of  time  at  a  fixed  pixel  location 
(236,211)  during  the  calibration  phase  along  with  quadratic  fit.  Data 
shown  is  typical  of  all  locations  throughout  the  field. 
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Table  1 :  Strain  data  obtained  using  fully  corrected  SEM  images  acquired  at  seven 
times  during  the  uniaxial  tension  loading  process. 
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Fig.  9:  Uniaxial  stress  versus  axial  strain  and  transverse  strain  for  uncorrected, 
drift  corrected  and  fully  corrected  image  correlation  data. 
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